Bogoliubov quasiparticle spectra of the effective rf-wave model 
for cuprate superconductivity 
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An exact-diagonalization technique on finite-size clusters is used to study the ground state and 
excitation spectra of the two-dimensional effective fermion model, a fictions model of hole quasipar- 
ticles derived from numerical studies of the two-dimensional t—J model at low doping. We show 
that there is actually a reasonable range of parameter values where the d^2_j,2-wave pairing of holes 
occurs and the low-lying excitation can be described by the picture of Bogoliubov quasiparticles 
in the BCS pairing theory. The gap parameter of a size Ad~0.13|1^1 (where V is the attractive 
interaction between holes) is estimated at low doping levels. The paired state gives way to the state 
of clustering of holes for some stronger attractions. 
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I. INTRODUCTION 

High-temperature superconductors are characteristic 
of the very short coherence length and small carrier num- 
ber, and may be in an intermediate regime between two 
limits, a BCS superconductor and a condensate of pre- 
formed bosons, and this observation may provide a pos- 
sible way of explaining the so-called pseudo-gap (or spin- 
gap) anomaly commonly observed in the lightly doped 
regime of cuprate superconductors ||]. Analyses of the 
two-dimensional (2D) t^J model at low doping have sug- 
gested a possible picture [||-||; under strong antiferro- 
magnetic spin correlations, low-energy properties may be 
described by the coherent drift-motion of the spin-bag- 
like quasiparticles where the rapid incoherent oscilla- 
tion of the 'bare hole' is eliminated as excitations inside 
the bag ||]. This picture immediately suggests an effec- 
tive fermion model, i.e., a fictions system of interacting 
fermions representing the hole quasiparticles ■ Since 
the attractive interaction acting between hole quasipar- 
ticles (or binding energy) is fairly strong ^0.8J in com- 
parison with the quasiparticle bandwidth ^2 J [||J|], one 
may hope that the effective model simulating this situa- 
tion would provide a possible explanation of some of the 
anomalous quasiparticle properties of cuprate supercon- 
ductors. 

The purpose of this paper is then to present an anal- 
ysis of the ground states and some excitation spectra of 
finite-size clusters of the effective fermion model at T=0 
K. We show that there is actually a reasonable range of 



parameter values where the c?-wave pairing of holes oc- 
curs and the low-lying excitation can be described by the 
picture of Bogoliubov quasiparticles in the BCS pairing 
theory [p^ . Also we show that the d-wave gap param- 
eter is nearly in proportion to but much smaller than 
the strength of the attractive interaction between holes. 
However, clustering of holes occurs for some stronger at- 
tractions, unlike, e.g., in the negative-C/ Hubbard model 
where the coherence length continuously decreases with 
the attraction strength and a condensate of tightly-bound 
real-space pairs becomes a good description 11 1. 



The effective fermion model considered here may be 
defined by the Hamiltonian 
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with a negative value of the density interaction V act- 
ing on holes in the nearest-neighbor sites <ij> 
n.i is the number operator at site i, and cj^^ creates a 
hole ('spinless' fermion) with momentum k in the cr (=|, 
i) sublattice; we have the two-sublattice model where 
the holes with up-spin (down-spin) are always in the |- 
sublattice (j-sublattice) and the numbers of up-spin and 
down-spin holes are the same. The single-hole dispersion 
in the 2D antiferromagnet is taken as the noninteracting 
band structure 

£k — 4tii cos kx cos ky -\- 2t2o(cos 2kx + cos 2ky) 

where in and ^20 are the second- and third-neighbor hop- 
ping integrals with positive sign. We employ the Lanczos 
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exact-diagonalization technique on finite-size clusters of 
up to 64 sites with periodic boundary condition; we can 
take clusters much larger in size than those feasible for 
the t—J model. Hereafter we refere to the momentum 
defined with respect to the 'nonmagnetic' Brillouin zone 
(see Fig. 1). 

II. GROUND STATES 

The ground states with zero total momentum are cal- 
culated as a function of the parameter values ^20/^11 and 
V/tii for all the clusters with 32, 36, 50, and 64 sites. 
We find that, irrespective of the size of the clusters, the 
ground states have the point-group symmetries in the pa- 
rameter space, as those shown in Fig. 2; there are some 
quantitative differences in the phase boundary among 
clusters but the basic features are the same. 

First, we note that two holes in the empty lattice 
form a dj;2_y2-wave bound pair for ^20/^11^0. 6— 0.8; this 
is simply because the two holes can gain kinetic en- 
ergy only by changing the spatial sign of the pair wave- 
function since a hole goes round the other with the pos- 
itive hopping parameter tu For t2o/^ii^0.6— 0.8, 
this mechanism of d-wave pairing works less favorably, 
and also in the k-space, holes tend to accumulate around 
k=(±7r/2, ±7r/2) where there is a node of the d-wave pair 
wave-function, so that the p-wave pairing occurs. There 
are two regions in the four-hole case; at ^20/^11^0. 6 and 
1^1/^11^4 we have the region where there are two d-wave 
pairs (as will be shown in Sec. IV). In the six- hole case we 
again find a parameter region where three d-wave pairs 
present. Clustering of holes occurs however for larger at- 
traction strength, suggesting the model to be phase sep- 
arated. This is evident in the calculated density correla- 
tion function {niUj) where (...) denotes the ground-state 
expectation value; we find that with increasing |T^|/iii 
the probability of finding a hole in the nearest and next- 
nearest neighbors of a site with a hole increases sud- 
denly across the level-crossing line in Fig. 2 and satu- 
rates rapidly to a constant value, and it also increases 
very rapidly around y/tii~— 2 even for <2o/^ii^0.5. 

We should note here that in the effective model the hole 
can jump over the other because of the hopping term to 
the second and third neighbors. This process however 
is forbidden in the t—J model. The stronger tendency 
to hole clustering in the effective model may therefore 
be suppressed by somehow implementing a 'conditional 
hopping' of the t—J model so as to gain in kinetic energy 
of holes. Such a correction however will be discussed else- 
where. Hereafter, we present our analysis of the model 
Eq. (1) in the parameter region where the clustering of 
holes does not occur (i.e., smaller |F|/iii values) and 
two holes in the empty lattice form a d-wave bound state 
(i.e., ^20/^11^0. 5). We thereby examine what are the 
low-energy excitations in this parameter region. 



III. SINGLE-PARTICLE SPECTRA 

We first calculate the single-particle spectral function 
defined as A(k, oj) = A"(k, -lo) + A+(k, cj) with 

A-{k,u;) = ^ |(V'^-i|ck.|Vo^)P'5(c. - Ei:'-^+E^) 

A+(k,c.) = ^ |(v^.^+^i4jOP^(- - e:^+' + o 

(2) 

which simulates the angle-resolved photoemission (PES) 
and inverse photoemission (IPES) spectroscopy. E^^ and 
are the ly-th excited eigenenergy and eigenstate {1^=0 
denotes the ground state) of the TV-hole system of the 
cluster, respectively. We add a small imaginary number 
ir] to oj and give a Lorentzian broadening to the spec- 
tra which otherwise consist of a set of 6 functions; the 
value 77/tii=0.05 is used. This spectra can in principle 
be compared to recent experimental data 

The calculated results for ^(k, ijj) are given in Fig. 3. 
When value is small the spectra simply refiect 

the noninteracting band structure, but for larger values 
of we find that the spectra are sharp near the 

Fermi energy but become broad away from the Fermi 
energy and momentum, reminiscent of the 'dressing' of 
the particle of correlated Fermi-liquid systems; the Bloch 
electron added into the correlated Fermi sea may be scat- 
tered out of its single electron level, leaving the system 
in the manifold of excited states. 

The calculated spectra also show how the attractive in- 
teraction affects the single-particle excitation spectrum. 
The gap-like structures can be seen in the spectra at (and 
around) kp, and the size of the gap increases with in- 
creasing attraction strength. The spectral feature, i.e., 
spectral- weight transfer appearing in both PES and IPES 
sides near kp and vanishing far away from kp, is noticed, 
which is consistent with the expectation of the BCS the- 
ory. This feature is absent along the diagonal {kx=ky) 
of the Brillouin zone and is consistent with the d^^-y^- 
wave pairing. The spectral-weight transfer also indicates 
a tendency toward smearing of the jump at kp in the 
momentum distribution function. 



IV. BOGOLIUBOV QUASIPARTICLES 

Now let us examine the validity of the Bogoliubov- 
quasiparticle picture for the low-lying states of the ef- 
fective model. We use the technique proposed in Re£ 
|]l2| , i.e., an exact calculation of Bogoliubov quasipar- 
ticle spectrum on small clusters, and see if this picture 
works in this model. We define the one-particle anoma- 
lous Green's function as 

G(k,z) = «^'l4T ,,g\j^ ^UM^) (3) 
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with the spectral function 

F(k, uj) ^ --ImG(k, uj + iri) 
n 

and its frequency integral 

where Eq is chosen as the average of and E^^"^. 
The hypothesis (see Ref. for details) that low-lying 
states of the clusters can be described by the microcanon- 
ical version of the BCS pairing theory then predicts that 
F(k, w)=Fk(5(w - Ey,) with Fk=Ak/2£:k, where E^, and 
Ak are the quasiparticle energy and gap function, respec- 
tively. Similarly, we define the two-particle anomalous 
Green's function as 

G(k,k',z) = (V'^lcLk.xC^ki ^ pjv 4'TgkTlV'^) 

Z ~ £1 -\- HiQ 

_n^kxnkT,^^, (4) 

z 

with nk(T=(V'o'|cj^^Ckcr|V'o')- Only the iV-particle sub- 
space is involved here unlike in Eq. (^). We define 
the spectral function _F(k, k',w) and its frequency in- 
tegral i^kk' as above. The hypothesis then predicts 
FCk^h' ,uj)—FkF^,d{uj—Ek—Ek'). By examining these 
two anomalous Green's functions, we can see if the low- 
energy excitations of the model are described by the BCS 
pairing theory. 

The spectral function F(k, w) is calculated for a 64-site 
cluster of the effective model and is shown in Fig. 4 for a 
number of the attraction strength. Noticing that the cal- 
culation is made by the subtraction of two single-particle 
excitation spectra [|l2|, we should first of all stress that 
nearly all incoherent spectral features in A(k, to) are sub- 
tracted out, leaving only the low-energy features (stem- 
ming from weakly-interacting Bogoliubov quasiparticles 
as shown below). 

We then find the following, all of which are consis- 
tent with expectations of the BCS pairing theory: A 
pronounced low-energy peak appears at kp and smaller 
peaks appear at higher energies for other momenta; the 
weights of the peaks are consistent with the BCS form 
of the condensation ampitude fk with a maximum at 
kp (see Fig. 5). The momentum dependence of F(k, w), 
i.e., the change in sign under rotation by 7r/2 and van- 
ishing weight along the kx=ky line, is a clear indica- 
tion of d^2_y2-wa,ve pairing. The size of the energy gap, 
which may be estimated directly from the positions of 
the peaks, increases with increasing value. With 

decreasing gap size, the peaks at momenta other than 
kp lose their weight as expected from the BCS theory. 
We note that, when the attraction strength V/tn ex- 
ceeds some critical value, these features are collapsed 
and the spectra completely loose the significance of Bo- 
goliubov quasiparticles, indicating that the clustering of 



holes takes place and the system seems to phase sepa- 
rate. The coherence length ^ estimated from the mo- 
mentum dependence of Fk decreases with the attrac- 
tion strength but still is more than twice as large as 
the nearest-neighbor lattice spacing even near the crit- 
ical V/tii value. 

We further compare the calculated spectra i^(k, a;) 
with the BCS predictions where we assume the quasipar- 
ticle energy £'k=\/ (ck — fJ-)^ + A^ with the gap function 
Ak=Ac((cos fca; — cos/cy) and choose the chemical poten- 
tial /I to guarantee vanishing ek~M af ^f- The value of 
A^; is then evaluated by fitting the positions of the low- 
energy peaks in F(k, w). We find that the fitted quasi- 
particle spectra are in a fair agreement with the exact 
spectra, which demonstrates the validity of the BCS pair- 
ing theory for low-lying excitations in the effective model. 
The estimated values of the gap parameter A^ are shown 
in Fig. ^. We find that A^; is insensitive to the doping lev- 
els examined and has the value A^— 0.f3|V^| until reach- 
ing the region of phase separation, which is much smaller 
than the bandwidth for Ek- This value also corresponds 
to Ad~130K if we assume \V\^0.8J and J~1300K [§. 

A consequence of (i^2_j,2-wave pairing may be seen in 
the point-group symmetry of the ground states, i.e., an 
alternation of the symmetry between Bi and Ai for the 
two-, four-, and six- hole states in a parameter region 
(see Fig. 2). This alternation, present also in the t—J 
clusters |g2|, suggests the picture that holes are added 
in pairs with (ia;2_j,2-wave symmetry. G{k., z) in Eq. (|^) 
(and thus F{]i, uj)) reflects the pairing symmetry via the 
point-group symmetries of ■0;^' and ipo^'^- Note that 
F(k, k',w) is defined entirely within the A^-particle sub- 
space and thus is not affected by this alternation; never- 
theless it shows the same indication of d^2_y2-wave pair- 
ing as F(k,uj). The calculated results for F(k, k',a;) 
are shown in Fig. |^. We again find that the size of 
the excitation gap increases with increasing and 
the k-dependence of the spectra clearly indicates dx2-y2- 
wave pairing. There are sharp peaks at low energies and 
broadened features at higher energies. As is the case 
for F(k, w), these high-energy features lose their weight 
rapidly with decreasing value, whereas the peaks 

at kp become sharp but remain finite with decreasing 

1/^11 value. These results are consistent with the no- 
tion of weakly-interacting Bogoliubov quasiparticles for 
low-lying excitations in the BCS superconductors. When 
the attraction strength V/tn exceeds some critical value, 
the spectra F(k, k',a;) again completely loose the signif- 
icance of Bogoliubov quasiparticles, indicating that the 
clustering of holes takes place. 

V. SUMMARY 

By using an exact-diagonalization technique on finite- 
size clusters, we have studied the ground states and ex- 
citation spectra of the 2D effective fermion model, i.e., 
a system of spin-bag-like quasiparticles corresponding to 
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doping holes, which is derived from numerical studies 
of the 2D t—J model at low doping. We have exam- 
ined the point-group symmetries and density correlations 
of the ground states of various size clusters with two, 
four, and six holes over a wide range of parameter val- 
ues, and have found that the results arc insensitive to 
the size of the clusters and can be summarized as the 
schematic phase diagram reflecting the pairing symme- 
try of holes. We then have calculated the single-particle 
spectra of the model in a parameter region and have 
detected some indications of the d^2_y2-wsMe hole pair- 
ing in the obtained spectral features. We have further 
calculated the one-particle and two-particle anomalous 
Green's functions and found that there is actually a rea- 
sonable range of parameter values where the dx2-y2-wave 
pairing of holes occurs and the low- lying excitation can be 
described by the picture of Bogoliubov quasiparticles in 
the BCS pairing theory. The gap parameter is estimated 
to be Ad— 0.13|y| at low doping levels. It seems however 
that with increasing the attraction strength the state of 
clustering of holes overwhelms this dx2-y2-wave pairing 
state before the picture of a condensate of tightly-bound 
pairs with the dx2-y2-waye internal structure becomes 
appropriate. 
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FIG. 1. Available k-points in the Brillouin zone for (a) 32-, 
(b) 36-, (c) 50-, and (d) 64-site clusters with periodic bound- 
ary condition. The dotted lines indicate the Brillouin zone 
defined for the two-sublattice model. 

FIG. 2. Schematic representation of the point-group sym- 
metry of the ground state at k=(0,0) obtained for the 32-, 
36-, 50-, and 64-site clusters with (a) two holes, (b) four holes, 
and (c) six holes; Bi {d^2_y2), E {p^ andpj,), and Ai (s) sym- 
metries appear in the parameter space. 



FIG. 3. Single-particle spectra ^(k,a;) calculated for vari- 
ous attraction strengths {V/t=—l, —2, and —3 for left, center, 
and right columns, respectively) and hole numbers (two, four, 
and six holes), ui is the energy of the hole. The spectra for the 
64-site cluster are arranged from top to bottom in each panel 
as k=(7r/2, V2), (^/4, ^4), (0,0), (7r/4,0), {n/2), (37r/4,0), 
(tt, 0), (37r/4, 7r/4), and (7r/2, 7r/4). Those for the 32-site clus- 
ter are as k=(7r/2, 7r/2), (7r/4, 7r/4), (0, 0), {■k/2, 0), (tt, 0), and 
(37r/4,7r/4). 

FIG. 4. Bogoliubov quasiparticle spectra F(k,cj) for a 
number of attraction strengths. The upper four panels show 
the results from the anomalous Green's function between 
two- and four-hole ground states for the 64-site cluster with 
t2oAii=0.3, and the lower four panels show the results be- 
tween four- and six-hole ground states for the 32-site cluster 
with t20 7^11=0.3. The arrangements of the spectra in each 
panel follow those in Fig. 3. 

FIG. 5. Attraction strength dependence of the condensa- 
tion amplitude calculated for the 64-site cluster with the 
filling between two and four holes (left column) and for the 
32-site cluster between four and six holes (right column). 

FIG. 6. Gap parameter Ad/tn as a function of the attrac- 
tion strength estimated from the calculated F'(k, w) 
for two-, four-, and six-hole systems. 

FIG. 7. Bogoliubov quasiparticle spectra _F(k, k',a;) calcu- 
lated for the 32-site cluster with six holes. The spectra for var- 
ious k points with k'=(37r/4, — 7r/4) are shown. t2o/tii=0.3 
and V/tii=—2.0 are used. 
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